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SUMMARY 


A computation package containing all equations and procedures needed in 
designing a high-power multicavity klystron amplifier has been developed. 

The rigorously derived three-dimensional relativistic axi symmetric equations 
of motion are used to compute the bunched current and the induced RF gap 
voltage for all interaction cavities except the input and second cavities, 
where the linear space-charge wave theory data are employed in order to 
reduce the computation time. Both distance-step and time-step integration 
methods are used to compute the Fourier coefficients of both the beam cur- 
rent and induced current. 


INTRODUCTION 

This report describes in detail a design computation package of a mag- 
netically focused high-power multicavity klystron amplifier. The theoreti- 
cal foundation is based on the investigation of field-electron interactions 
of a relativistic klystron described in a previous report (ref. 1). The 
computation method is based on the large-signal disk-electron model which 
can be programmed and run on a large-capacity high-speed computer. The un- 
modulated electron beam which enters the input cavity is divided into a 
large number of deformable disks of electrons, and each disk is subdivided 
into an equal number of rings of electrons which serve as elementary charges 
in the mathematical model. Thirty-two disks per one space-charge wavelength 
and three rings per disk were used in the computation model. Each of the 
ninety-six rings is followed through from the input cavity gap, to the 
intermediate interaction gaps, and to the drift-tube spaces until the output 
gap is reached. The induced current and the RF gap voltage of the various 
cavities are computed on the basis of the theory of Shockley and Ramo. The 
fundamental bunched current and the induced RF gap voltage of the various 
cavities are computed, first, by performing the distance-step integration of 
the equations of motion to yield the precise data of the position and phase 
of the individual rings and second, by evaluating the Fourier coefficients 
of the bunched beam current and the induced RF gap voltage by using the 
time-step integration method. Beam loading and cavity detuning effects are 
appropriately dealt with. Space-charge force Green functions, Sfr 
and Sfp, are modified to take into consideration relativistic velocity 
of the beam motion. Since the RF gap voltage and the beam loading parame- 
ters, Gh s ^b s> S'"® interrelated, an iteration process is em- 

ployed t6 determine’the voltage modulation index for the various gaps. 
However, in order to reduce the computation time the linear space-charge 
wave theory data were used to design the input and second cavities, where 
small voltage modulation index values will be used. 


Material, not previously published, presented at the International Electron 
Devices Meeting sponsored by the Institute of Electrical and Electronic 
Engineers, Washington, D.C., October 11, 1971, and the Solar Power Space 
System Workshop sponsored by NASA Lyndon B. Johnson Space Center, Houston, 
Texas, January 15-18, 1979. 


This report consists of three main parts: 

(1) All input data and initial computation formulas needed to start the 
computer program 

(2) All working equations used in various stages of computation 

(3) Design procedures described in detail, including the method of 
iteration in the computation of the voltage modulation index of the 
various gaps 

In addition, four appendixes are included at the end of this report so 
that the origin of some of these working equations can be traced and checked. 

INPUT DATA 


Design Parameters 

The following are the parameters and data to be specified by the 
designer: 

(1) Klystron operational frequency, fQ, in hertz. (This is also the 
input drive frequency of the klystron.) 

(2) Input drive power, P-jp, in watts. (Matched condition is assumed.) 

(3) Direct-current beam current, Iq, in amperes, which is related to 

the relativistic perveance by the equation 
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where Kg = Io/(Vq)^^^ is the nonrelativistic perveance, and 
V^q = mgc^/lel = B.llxlO^ volts is the equivalent beam voltage. 

(4) Direct-current beam voltage, Vg, in volts. 

(5) External applied dc magnetic field as measured along the z-axis, 

Bg, in weber per meter squared. (I Wb/m^ = 10 000 G). 

( 6 ) Tunnel diameter, 2a, in meters. 

(7) Direct-current beam diameter at the entrance to the input cavity, 2b, 
in meters. 

( 8 ) Diameter of the cathode, in meters. 

(9) Physical gap length of the s^^ cavity, 2 JI 5 , in meters, with 


s = 1 , 2 , 

2)1/1, and 2 ic. 


that is, for a five-cavity klystron, 2ii, 2 «- 2 » 2 ) 13 , 


(10) Drift-tube length, Lg, in meters, as measured between the 
centers of the s^^ and the (s ± 1 )'^*^ cavities; that is, 

Li, L 2 , L 3 , L 4 , and L 5 . It is to be noted that the 
drift-tube length may be determined initially by the small-signal 
space-charge wave theory data given by 


'■s 


sp,(s) 


= ( 2 m + 1 ) ^ 
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'po 


where m = 0 , 1 , ... and 
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electronic wavelength, y ~ 
space-charge wavelength 

5 

relativistic dc beam velocity, 5.93x10 m/sec 


plasma radian frequency. 




relativistic velocity reduction factor. 



nonrelativistic dc beam velocity. 


5.93x10^ ’’’/sec 


(11) 00 = 2irfo/uo- 

(12) k = 2irfo/c. 

(13) Nonrelativistic cyclotron radian frequency, wqo = 

(14) Unloaded Q of the s^*^ cavity, Qu^s? ^I^^t is, Qu,l. Qu,2» 

Qu,3’ Qu,4’ Qu,5-.. 

(15) External Q of the s^^ cavity, Qext,sJ that is, Qext,l ^nd 

Qpyi" R (Qpyh Opxi" 3nd Qpy 1 “ 4- ~ • 

(16) Characteristic impedance of the cavity, (R/Q)^; that is,* 

(R/Q)l, (R/Q) 2 , (R/Q) 3 , (R/Q)4, and {R/Q)5- 

(17) Frequency-tuning parameter of the s^h cavity 5u s = 

(ttuned * t)^esonant) /t^espnant^ that is, «u, 2 ’ '^u,3» 

^u,4 ' u , 1 3nd ^u,5 ~ , 

(18) Field-shape parameter of the stn cavity, Hg; that is. Hi, H 2 , 
H 3 , H 4 , and H 5 . 

(19) Number of disks in one space-charge wavelength used in the design 
model (here 32), N, and number of rings in each disk (here 3), R. 


Physical Constants Used in Computation 

The following physical constants (expressed in SI units) are needed in 
the computation: 

(1) Velocity of light, c = 3x10^ m/sec 

(2) Dielectric permittivity of the vacuum, gq = (1/36 ir)xlO“" = 
8.854x10-12 F/m 

(3) Magnetic permeability of the vacuum, pQ =_4 ttx 10-' H/m 

(4) Rest mass of the electron, mo = 9.108x10-21 kg 

(5) Magnitude of the electronic charge, |e| = 1.602x10-1" C 

( 6 ) Ratio of electron charge to rest mass, le/mgl = 1.759xloH C kg-1 

(7) Equivalent beam voltage, Vgq = moc2/|e| = E.llxlO^ V 
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WORKING EQUATIONS 
Normalized Equations of Motion 


The normalized axial equation of motion is 


= 




1 ]' ^ 




NR \ 0) 



3ea 


[1 - (ka) C ] 





r .^2 

e> 




( 1 ) 


The normalized radial equation of motion is 


= 


H 1 

“s ®s ^0 sinh[H^a(jig7aT] 


4 — (3 a)' 
a ^ e ^ 


[1 - (ka)2 p2] (pP 






The independent variables are 5, P0,eq,j’ defined as 

C normalized axial coordinate, z/a 

Pn,eq,j normalized equivalent charge center of ring charge 

element at entrance to input cavity 
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-j entry phase of disk charge element where j = 1, 2, 

’ 32 (=NR) in a 32-disk 3-ring model, 2irj/N 

» • 

The independent variables are p, p, and c, defined as 

p normalized radial coordinate, r/a 

p(5,Pn ea i’*n i) normalized radial velocity component of ring 

charge element, where 0 = ut, dp/de 

i(C,Pn eo i’'l’n i) normalized axial velocity component of ring 

charge element, d5/de = z/toa 

4>(€,pQ^eq^ j) phase variable which defines phase position of 

ring charge element relative to RF gap filed as a 
function of displacement and initial state of enter- 
ing electron beam, e - (ee^)? 

The following are the relationships or equations needed to compute equations 

(1) and (2): 

(1) Relativity parameter: 



41a tunnel flux, ira^ Bq 

4'c cathode flux, ur^ Bq 


r^ cathode radius 

Bp(p,e), B^(p, 5) radial and axial components of dc magnetic fields 

given by eqs. (3) and (4) in the section Calcula- 
tion of External Direct-Current Magnetic Fields 
Bp and B^ 
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(2) Four field-shape functions, F^, Fp, G^, and Gp: 
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Gp(5,p) 


■I 


Jl(px„) sinh(p„ H^a) -ji 
Pn * 


n.1 


^ sinh(p^ - H^a) / 
Pp - 


-PpICl «s 

e for Id > — 

d 


(3) Two Green function space-charge forces, and 





where Jq and are Bessel functions of the first kind, Xp is the n^*^ 

root of the Bessel function. 


sign 5 - ^0 " 


r 1 for 5 > Cq 
\- l for C < Cq 


Pp- /’•'i - 
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The notation (5 - Cf,) implies that the value inside the parentheses must 

R 

be computed by using the retarded positions defined by 


d = 




r ) 
eq,a' 




or 


(«'=) = »eq,a>' 

(see fig. 1); that is. if we wish to calculate the force acting on the 
charge located at c at time t (i.e., e) due to a charge at we 
must take the value of charge at time t - |z - zo|/c, or more precisely, 
t - d/c, where d is given above and c is the velocity of light. 

Calculation of External Direct-Current Magnetic Fields Bp and 

With reference to figure 2, we note that the radial and axial compo- 
nents of the external dc magnetic field produced by an electromagnet are 
given by (appendix A) 



Figure 1. - Computation of space charge forces between 
disks (rings). 
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Bp(C,p) 


r '-sc 

/ Kj^(ka^)Ij^(kap)sin(kac)sin 

*^r\ 


= B, 


n=l 


/ Kj^(ka^)sin ka 


(3) 


dk 


B^(5,p) 


^ y Kj(ka, 

•/n 


)lQ(kap)sin(kac)sin ka dk 


= B. 


n=l 


^ J K^(ka.)si 
*'o 


(4) 


sin ka dk 


where Iq, Ii, and are modified Bessel functions of the first 
kind; Bq is the magnetic field along the z-axis; and a^ is the mean 
radius of the layer. It has the following simple relation with the 

diameter d of the wire used: 


3i = ai + (i - l)d 


where a^ is the mean radius of the first layer of the solenoid. 

Calculation of Fundamental Bunched Current for Third and Succeeding 
Cavities, s = 3, 4, ..., Including Output Cavity 

The fundamental beam current when normalized to the dc beam current 
Iq is given by (appendix B) 


+ B? e ^ 


(5) 


s s 


where 


NR=96 


^s “ NR S *j(PO,eq’''’0’^s^ 


( 6 ) 


j=l 


NR=96 

^s=m Z ^,5 

j=l 


(7) 
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Figure 3. - Identification of various interaction gaps and drift-tube spaces in 
five-cavity klystron. 


A. ^ = arctan — (8) 

b.s 

where the bar over the symbols Ag and Bg indicates they are normal- 
ized Fourier coefficients, and 

ri for ? = (i.e., z = z^) 

^’^s \o for 5 it (i.e., z 4 z^) 

in which Cc (i.e., Zg) is referred to the midplane of the s^^l^ 
cavity gap (fig* 3) with 



for i = 1, 2, 3, . . . (9) 


For example, for the fifth cavity, s = 5 (see fig. 3), 





a 


The equivalent charge center pg eq,i (i.e., ro,eg,j) 9^ 

charge element is determined by the design data sn6wn in figures 4 and 5 
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-r*/' RING OF 
I n'^DISK 



^0 

*0 


(a) UNMODULATED DISKS BEFORE 
ENTERING INPUT CAVITY GAP. 



‘s 


(b)TWO MODULATED DISKS AT 
z - Zj AND t - tj WITH POSSIBLE 
SUBSTANTIAL DEFORMATION. 


Figure 4. - Unmodulated and modulated beam disks. 


rg_ 3 = bR 





Figure 6. - Identification of ring charge eiement of three-ring 32-disk model. 


11 




5 ^ O.eqg _ b R 

0,eq,l a a ^ 


_ ^0,eq,2 _ b _R 
PQ,eq,2 a a 




y 


( 10 ) 


P0,eq,3 


0,eq,3 

a 


b R 


VI 


J 


In order to facilitate the computation, we introduce the equivalent charge 
center with the general index j and the index r to designate the ring 
of a given disk. Then, for a three-ring model (R = 3), r = 1, 2, and 3. 
Furthermore, if we introduce a second index d to designate the number of 
disks under consideration, with d = 1, 2, 32 for a 32-disk model, we 

can replace the index j in Po,eq,j by 


j = 3(d- 1) + r 


for 


r = 1, 2, and 3{R) 
d = 1, 2, 32(N) 

For example, with reference to figure 6, we note that the 95l-b ring 
charge element (j = 95) is identified as the second ring (r = 2) of the 
32'^° disk (d = 32); hence 


j = 3(32 - 1) + 2 = 95 


On the other hand, since each disk and each ring become distorted and de- 
formed after the beam has passed through the input modulation gap, in order 
to identify and follow the movement of a given ring charge element the 
following equation is used instead: 


'"s,eq,j '"0,eq,j 



dt 


or 


*^s,eq,j 


'~s,eq,j 

a 


^0,eq,j 


+ 



ap de 


( 11 ) 
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r 


Calculation of Beam Parameters Sb,s ^b,s 

The beam conductance beam susceptance Bt>,c needed 

to compute the gap voltage modulation index cxg are determined by 
(appendix B) 






cos Xg 


( 12 ) 





sin 


(13) 


where 


I 


0 


^0 ^ T" 

U Vq 


B 

c c 

x^ = arctan(2QT ^6-r ^) + arctan arctan — 

s 1 ,s T,s j 


(14) 


with Ag and Bg given by equations (6) and (7), Cg and Dg given by 
equations (16) in the next section, and ag given by equations (18) in 
the section following that. 

Calculation of Induced Current I-j^g 

The normalized fundamental induced current is given by (appendix C) 

i ^ s _ + q 2 g '^'^i,s (15) 


where Cg and Dg 


are the normalized Fourier coefficients given by 



-2ir Ha 


MNR sinh 



M=36 


I z 


m=l 


all k 


NR=96 



j=l 




‘•’m^PQ.eq, j’*0'^eq,j’^k’*s^\ 
sin 't*m(Po,eq, j’^O’^eq, j’^k’^s V 


(16) 


13 


II 


and 


A . = arctan — 

1,S jr 


(17) 


with 




_ fl ^ 

k 1 0 for 5 ^ C. 


♦ 2) i 5 k i 's " " 2) 


It is to be noted that (1) the two equivalent charge centers PQ,eq,j 
and p pq j are determined in accordance with rules detailed in’the sec- 
tion CaTcuiation of Fundamental Bunched Current for Third and Succeeding 
Cavities, s = 3, 4, Including Output Cavity and (2) the two normalized 

velocity components pj and Cj are to be computed as functions of 



^ ' ^^0,eq,j’*0’Pk,eq,j’^k’*s^ 


Determination of Voltage Modulation Index ag 

The gap voltage modulation index for the s't*^ cavity is given by 
(appendix D) 


“s 






2 —2 

^ + D 
s s 


+ 4Q 


T,s'^T,s 


and the RF gap voltage is given by 


V 

g,s 


''o“s 



where 


G 



0 

0 


(18a) 


(18b) 
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1 


( 19 ) 


T,S - 1 ^ 1 ^ r 

^u,s ^ext,s 


6 - 6 + 

^T,s - ^u,s 2 \Q 


( 20 ) 


IT, 


e = -arctan(2Q-p 6-p ) - arctan 

S I 9 S I ) S p 


( 21 ) 


and TTs and Us are given by equations (16); Gb,s and 8^ c are 
given by equations (12) and (13); and Qu,s> Qext’s’ ^nd (R/(3)s are 
design parameters. 

COMPUTATION 
Input Cavity, s = 1 

Step 1 : Perform the following computations: 

(1) Compute the voltage modulation index ai by the relation 


“1 






Tr 


'0 


( 22 ) 


where P-j^ is the input power to the first cavity. Assuming that the 
input cavity presents a perfect match to the signal generator, we note that, 
since Q]oad = Qext,l’ 


^T,l 


+ G 


^u,l 


b,l VQ 


(22a) 


If the beam loading effect is neglected, Q-p^^ = Qu^i- 

(2) Compute the axial and radial beam coupling coefficients by using 
small-signal space-charge theory data given in the following: 


sin e 


M 


z,l 


e 




g,i 


(23) 


15 


where e,q^ 2 ^ is the gap half-transit angle defined by 


1 = B„a — 
g,l e a 


V,1 


/l^(Yb) - I^(Yb) 

lo(^a) 


( 24 ) 


in which Iq and are the modified Bessel functions of the first 
kind, and 


y = 


/i 


- 


with 


R. = 




(25) 


Step 2 : Use the voltage modulation index obtained in step 1 to 

set up the velocity modulation of the beam, and by integrating the equa- 
tions of motion (eqs. ( 1 ) and ( 2 )) compute the bunched current at the second 
cavity gap. The initial conditions for the electrons are 


z . Zg . 0 


for t = t. 


r = r, 


0 ,eq, j 


or 


5 = 0 


for 0 = e 


P = p 


0 ,eq,j 


0 


and 


"z = "O 


u = 0 
r 


for t = t. 
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or 


5 = 




r for 


p = 0 


e 


0 


Second Cavity, s = 2 


Step 1: Compute the fundamental beam current Ih ? by using for- 

mulas derived from the linear space-charge wave theory! This is given by 
(ref. 2) 


I 



2J, 




Tf 



(26) 


where 


Bessel function of first kind 
given by eg. (22) 

Mz^l given by eg. (23) 

given by eg. (24) 
given by eg. (25) 

Ab,2 -(eea)(Li/a) 

1-1 length of drift tube between first and second cavities (see fig. 3) 

toq reduced plasma freguency, Fup 

o)p plasma freguency 

F space charge reduction factor (determined from relation between 

Bgb and a/b) 

(For cylindrical beams in cylindrical tunnels, data on F can be found in 
ref. 2, fig. 3, p. 106, for instance.) 

Step 2 : Compute the induced current I-j 2 and induced voltage 

Vg^2 and hence obtain the gap voltage modulation index 02 - 

(1) The induced current (fundamental) is computed by small-signal 
theory formulas: 


1,2 

0 


= Mz,2Mr,2 





(27) 
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where 



sin 0 


g,2 


g,2 


in which 



0 0 
'^e 



a 


By equation ( 25 ), can be 

obtained: 



^•,2 

'*2 

■ ^b,2 “ ®e® a 


( 2 ) The fundamental 
(appendix D); that is. 

induced gap voltage is given by equation 

(D6a) 

V 

c 


|^•,2PT,2 (q)^ j62 

0 

( 28 a) 

1,2 - - 

f/l + 

where 





= A- , 

, - arctan( 2 Q^ 2^j 2^ 



— A L /- 

b,^ 

, - arctan( 2 Q-|. 

( 28 b) 


^T,2 

1 

" "*>'2 W 2 

( 28 c) 


^T,2 = 


( 28 d) 

S,2 ' 

^b,2 

^0 

Gq 

— cos[arctan( 2 Q, p6j 9)] 
»'■ 



^b,2 

1 

( 28 e) 


^0 

’2 |/l , 4Q^_2*^_2 


18 



^b,2 


b,2 


— sin[arctan(2Q-r ,)] 


b,2 


^^T,2^T,2 


and 




G 

So =7^ 


Hence the gap voltage modulation index 02 is 


(28f) 


“2 — 7 


,2 


= G, 


i,2 


^T,2 H 2 


f/l + 4Q 


t,2^T,2 


(29) 


We notice that 02 is a function of the beam loading parameters 2 
and B |3 2 which are, in turn, a function of the gap voltage. Hence,’ to 
compute’ « 2 , we require an iterative process. This is described in detail 
in the following: 


(a) Initially, we find 02 by computing Qj 2 6j 2 by using 

data of Gjj 2 ^nd B| 3^2 derived from small-signal theory: 


’b,2 


sin 9^ / sin „ 

0 q,2 I g,2 

7 e' I e 


cos 6 


9,2 


9,2 


9,2; 


2 ”^z,2^^z,2 " %,2^ 


(30) 


G„ /cos 9„ „ 

, 0/ q,2 

^b,2 - ■ 7“ I "9 
’ \ 9,2 


<”z,2 - “g,2> 


(31) 


The value obtained from equation (29) is given the name a2(0)’ where the 
subscript (0) denotes the initial data for 02 . In this manipulation, the 
values of Qj^2 6j 2 obtained will be assigned the names Qj^2(0) 

and 6-]-^2(0)' 

(b) In the next step, the values of a2(0)> Gt, 2(0), 3>^d <ST,2(Q) 

used to compute G[j 2 and B^ 2 fnom the large-signa V equations’(z8e) 
and (28f), respectiC'ely. The 3ata so obtained are given the name G^ 2(1) 
and Bh ?c\\, where the subscript (1) denotes the first iteration datl.' 
These new data of 2(1) and B^ 2(,1) are used to compute Qj 2(1) 
and «j^2(l) eqiiations (28c) Sno (28d), respectively, and h^nCe, to 

obtain a new gap modulation index, called «2(1)* 


19 



(c) After the first iteration process, we obtain an RF gap voltage 
Vg,2{l) given by 


V 


g,2(i) 


'^0“2{1) 



(32a) 


where 


®2(1) = ‘b,2 - 

and this is to be compared with Vg WO)’ which is defined in accordance 
with data obtained from the initial ueration process (a); that is. 


%, 2 ( 0 ) 


'^ 0 “ 2 ( 0 ) 



(33a) 


where 


®2(0) " ^b,2 ®'"‘^^^'^[^^T,2(0)'^T,2(0)] 

In order to compare two complex voltages, we must compare their magnitudes 
and phases separately as follows: 


“2(1) " “2(0) I - ^1 


and 


|®2(1) ®2(0) I - ^2 

where and two small specified values. If both of these 

inequality conditions are satisfied, we go to design step 3; otherwise we 
repeat procedures (a) and (b). This time the iteration process is initiated 
by using a^(i) obtained in (b), and after the process has been repeated 
n times, finally the following inequality conditions are simultaneously 
satisfied: 


|“2(n) ~ “2(n-l)| ^ ^1 


and 


®2(n) ■ ®2(n-l)| ^ ^2 


Then, a2(n) is the finalized value for 02 * The word "simultaneously" 
requires an explanation. It simply means that both inequality conditions 
must be satified at the same time. 
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step 3: The finalized value of ap (=ap(n^ described previously 

is used to set up the velocity modulation for gap 2. Thus the distance-step 
integration of the equations of motion is continued through the second 
drift-tube space L 2 until the third cavity gap is reached. 

Third Cavity, s = 3 

Step 1: Compute the Fourier coefficients A'^, and C 3 , F 3 of 

the fundamental beam current Ij, 3 and the fundamental induced current 
I-j 3, respectively. The relevant; equations to be used are (5) to ( 8 ) for 
I|3’3 and (15) to (17) for I-j 3. It is to be noted that the two veloc- 
ity components, oj and 5 j, and the phase factors, <i)j and 4 >ni, are 
determined by the'^distance-step integration of the equations of motion ini- 
tiated in step 3 of the last section. 


Step 2: The Fourier coefficients A 3 , B 3 and C 3 , D 3 , obtained in 

step 1 are used to compute the gap voltage modulation index 03 by using 
equation (18a). However, we note that 03 , like ap, is a function of 
the beam loading parameters Gj^ 3 and Bj, 3, which are, in turn, a 
function of the RF gap voltage ’Vg 3. Therefore 03 can only be deter- 
mined by means of an iterative process very similar to that described in the 
last section. Initially, the small-signal beam loading parameters, given by 
equations (30) and (31), are used to compute Qt,3(0) 3(0) 

using equations (19) and (20), where 3 = “’ 3 a design 

parameter. It is to be noted that, when equations (30) and (31) are used to 
compute Gk 3(0) and Bb^3(o), the gap transit angle 6 g should be computed 
by using the naif-gap length next step, the^values of a 3 (o)» 

Qt, 3(0), "^T 3(0) obtained above together with the four Fourier coeTfi- 

cients obtained in step 1 are used to compute the beam loading parameters 
^b,3(l) 3 nd 85 ^ 3 ( 1 ), by using equations (12) and (13) in this instance; 
“3(1) computed from equations (18) to (20). The iteration process 
described in the last section is followed through until a finalized value 
o 3 (p^ is obtained. 

Step 3 : Here, again, the finalized 03 (= 03 /^)) value is used 

to set up the velocity modulation for gap 3, and the distance-step integra- 
tion of the equations of motion is continued through the third drift-tube 
space L 3 until the fourth cavity gap is reached. 

Fourth Cavity, s = 4 


Computations needed to design the fourth cavity follow exactly the same 
procedures described in the last section; however, great care must be exer- 
cised in choosing the appropriate design parameters for the various for- 
mulas. For instance, in computing the small-signal theory beam loading in 
the initial interation process, X-s should be used to compute the half- 
gap transit angle 0 q, etc. 

The design of tne fifth and sixth (if any) cavities follows in like 
manner. 


Output Cavity 

Step 1 : Set the index s in various formulas in cavity 2 to its 

appropriate value; that is, for a five-cavity klystron, set s equal to 5. 
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step 2 : Compute the RF gap voltage Vg^s from equation (18b), where 

the voltage modulation index for the output^cavity is obtained from 

an iterative process described in the sections Third Cavity, s = 3 and 
Fourth Cavity, s = 4. However, since the output cavity is coupled to an 
external matached load, some appropriate value (a design value) must be 
assigned to Qext,s when computing Qj^s equation (19). More- 

over, since the output cavity is tuned to the input driving frequency, 

6 ^ 1 , s can be set equal to zero when computing from equation ( 20 ). 

Step 3 : The following additional computations are needed to complete 

the design of the output cavity: 

(1) The fundamental conversion beam power, that is, the fundamental 
power taken from the beam by the output cavity, is given by 

P = Re I . , V* . 
conv \ 2 ) ind,s g,out 



cos[arctan(2Q^ ^)] 


This is maximum when cos[arctan( 2 Qj^ 56 y 5 )] is equal unity, that is, 
when the following condition is satisfied: 

arctan( 2 Qx ^ 6 -r ^) = 2 irk for k = 0 , 1 , 2 , ... 
i 9 S I ^ s 


(2) The conversion efficiency rtconv* electronic efficiency 

rig, is given by 


*^conv 


^e = 


conv 


(3) The actual power delivered to the external matched load Gl is 
given by 


P 


out 


( 2 ) '^g,s’load 

(0 


„2 2 

Vs 


^ext, s 
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(4) The power gain of the klystron is defined by 


• L/U O 

power gam = -r 

^in 

(5) If Pio§s,s dissipated in the output cavity wall, 

the output circuit’eff iciency is given by 

_ 1 _ power dissipat ed in cavity 
^cct “ power supplied by beam 

where power dissipated in the cavity is computed by 


loss,s 7R 


g,s 


= P 


^ext,s 


sh, s 


out Q 


u,s 


^sh s "'5 shunt resistance due to wall and output coupling losses 
alone. 

(6) If P^q] is the power recovered in the depressed collector. 


out 


out 


net - bean, power - 



APPENDIX A 


DERIVATION OF EQUATIONS FOR CALCULATION OF Bp(r,z) AND Bz(r,z) 

OF DIRECT-CURRENT MAGNETIC FIELD 

With reference to figure 2, we note that the vector potential A(r,z) 
due to a single-layer close-wound solenoid of radius a-j and length 
l_sc is given by (ref. 3) 


Ag(r,z) 


= C 



K 2 (xa^.)l 2 (^'")cos xz sin 



dx 

X 


where Ij and K 2 are modified Bessel functions of the first and 
second kind, respectively. The required magnetic field B is obtained by 
taking curl of A to obtain 


B(r,z) = V X A = a^ 


_-rli ('■%) =2 


7TF 


(Al) 


This gives 


sAq 

= - — 


= C 


^ K2(^3^)l2( 


»(¥) 


(xr)sin xz sinl^^— )x dx (A2) 


and 


^0 


. C 


f 


) lQ(xr)cos xz sin I -4^ 1 x dx 


>■{¥} 


(A3) 


In order to evaluate the constant C, we require that 


B^(0,0) = Bq for r = 0 and z = 0 

where Bq is the dc magnetic field along the z-axis, a design parameter. 
Thus, by letting r = 0 and z = 0, we have from equation (A3) 
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c = 


J Ki(>^ai)sin 


(A4) 


X dx 


When the value of C is substituted into equations (A2) and (A3), we get 


Bp(p;E) - 


f 

•/ n 


(xap)sin(xa?)si 


■(¥) 


X dx 


J K^(xa.)sin 


X dx 


(A5) 


and 


B^(p,0 = 


f ^i<“i 

•'A 



)lQ{xap)sin(xae)sin I X dx 


f Kj(xa.)sin 

•'o ^ 


X dx 


(A6) 


Equations (A5) and (A6) give us the radial and axial components of the dc 
magnetic field due to a close-wound single-layer solenoid. If the solenoid 
consists of N layers, the total magnetic field generated by this electro- 
magnet is obtained by summing all contributions of the N layers; that is. 


Bp(p,?) = 


/ Kj^(xa^.)I]^(Aap)sin(xa?)s1n x dx 


i=l 


J K^(xa.)sin 


X dx 


(A7) 


B^(p,0 = 


J Kj^(xa^)lQ(xap)sin(xa5)sin X dx 


j K^(xa^)si 

J r\ 





sm I I X dx 


(A8) 
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where Iq is the modified Bessel function of the first kind of the zero 
order, I]^ and are the modified Bessel function of the first and 
second kind, respectively, of the first order, and a^- is the radius of 
the layer of the solenoid, with i = 1, 2, ..., N. 
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APPENDIX B 


DERIVATION OF EQUATIONS FOR CALCULATION OF BUNCHED 
CURRENT AND BEAM LOADING 

The computation of the bunched current is based on the charge conserva- 
tion principle. If the beam is divided into a number of disks of electrons 
and each disk is further divided into elementary charge rings (figs. 6 and 
7), an elementary ring ir(r5 o ~ *^1 o)'^^0 containing a charge Aqg 

at time to at a later time t becomes ’aq in an element of the charge 
ring aS az. Thus, by the charge conservation principle, we obtain 

" Aq(r,z,t) 

or, more precisely, 

sz (B1 

Since the rings are deformed and are permitted to change their axial and 
radial dimensions as well as their shape according to the total force acting 
on them, the right side of (Bl) is employed to account for this fact; aSj 
is the surface element of the ring at time t. Equation (Bl) can be 

written as 


Pj(r,z,t) = pQ(rQ,Zp,tQ) 



AS. 


AZ 





u^(r,z,t)ASj A* 


(B2) 



Figure 7. - Computation of charge conservation in deformed disks 
(rings). 
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where az = At, a< 1) = w At, and aJzO = -pq'^zO- density of the 

current carried by the charge element aqj at time t and some distance 
z from ZQ is given by 


AJj(r,z,t) = Pj(r,z,t)Uj.(r,z,t) 


= Pj(r,z,t)u^(r,z,t)T^ + Pj(r,z,t)u^(r,z,t)a^ 
and the axial and radial components of the current density are given by 


AJ^j.(r,z,t) = Pj(r,z,t)u^j.(r,z,t) 


■ '^zO^'^0’^0’^0^ 


’('"o.o ''no). 




A4> 


(B3) 


and 


AJ^j(r,z,t) 


‘^z0^'"0’^0’^0^ 


’"Vo,o ^ 


2 \ 

1.0/j u^j(r,z. 


t) A«t), 




u^j(r,z,t) a4> 


(B4) 


respectively. 

Assuming that there is no net current flow in the radial direction, we 
obtain the total current passing through a given plane perpendicular to the 
beam at an arbitrary distance z and the time t by summing all contribu- 
tions from the charge elements that pass through this plane simultaneously 
at the time t. Hence 


j=NR 


j=l 


j=NR 




Z . / 2 _ ^2 \ 

®z,Zgy0,0 ■'.0/j A<|)(rQ,tQ;z^,t 


(B5) 


j=l 


where 6 is the Kronecker delta notation defined by 


5 =1 for z 


z 


s 
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r 


6 =0 for z 4= z 

and s is designated for the number of the cavity gap concerned. The sum- 
mation process requires further explanation. Since in a high-density modu- 
lated electron beam each disk of electrons is deformed and distorted in 
shape, only a portion of the charge may cross a given cross-sectional plane 
at a given time; and at a given time, charges belonging to different disks 
may cross a given cross-sectional plane simultaneously. Only the ring 
charges whose arrival time at z = coincides will contribute to the 
bunched beam current. 

In the next step, we expand the bunched current, equation (B5), in 
terms of Fourier series as follows: 


n=l 

where 



This can be evaluated by using equation (B5) for ib(zs»ts) to obtain 






E 


j=i 






(B6) 


and in like manner, 



1 

7T 







d4> 






2 \ /cos nV 

''4Vsinn. 


d*o(ro>^0’^o) 


(B7) 


Further evaluation of the coefficients Aq, Ap, and Bp requires some 
thought. Let us consider the simpler case of the evaluation of Aq first. 
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With reference to figure 6, we see that since by our choice each charge disk 
which enters the input cavity gap at equal time intervals is unmodulated, 
all the rings in a given disk will enter the input cavity gap simultaneously 
at a given time. Thus, if we use a double subscript to identify each of the 
NR charge rings instead of using a single index j, we can write equation 
(B6) as follows: 

NR 

^0 = " (?) '^z0^’"0’^0’^0^ 2 2 X (’"0,0,r " ’^i,0,r) ^*0 

all Hq n=l r=l 


where the subscript r designates the ring of the disk, and the subscript 
n designates the disk number. With this new convention, we can represent 
A<|)o by 2ir/N; hence 



(0'«(r) Z E ( 


n=l r=l 


0,0, r 



where 



N 




0 ’ ^ 0 ’ ^0 


) 


N 


Z 


b 


2 

n 



’"'^zO^'^O’^O’^O 


)(Nb^) 


= -I 


0 



(B8) 


is obviously the square of the beam radius b. This is the correct result 
for the coefficient Aq. The evaluation of the coefficients Ap and Bp 
is not as simple since we are dealing with multiple value functions in tne 
bunched beam. In this case, a single index representation may be least con- 
fusing. Thus, we shall represent a4>q by (2 tt/NR). Using this conven- 
tion, we may write equation (87) as 



2tt 

NR 


‘^zO^'^O’^O’^O^ 



(B9) 


where 4 >j = «i)j(rQ,tQ;r,Z 5 ,t 5 ) is a function of the radial position of the 
charge ring on crossing the constant Zg plane at time tg, and hence. 
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it is to be determined by its radial position at the entrance rg and its 
entering time tg. In computer programming, it is more convenient to 
identify an individual ring by its equivalent charge center rgq or 
rg instead of its two radii, the outer radius rn and the inner radius 
r^-. Using this convention, we write equation (69) as 


( 





6 


Z,Z 


s 



cos nct>j(rQ,tQ;r,r,Zg,tg) 
sin n<tj.(rp,tQ;r,Zg,tg) 


(BlOa) 


where Iq = ^be dc beam current. In the expression for the 

equivalent charge center t"o,eq,j» ^be subscript 0 is referred to the 
unmodulated beam at the entrance to the input cavity gap. It is to be noted 
that, if each ring of the disks is divided into R rings, then, in accord- 
ance with the double index notation described previously, we have the 
following relation: 


*"0,eq,r ~ 


= r, 


0,eq,r 


n=N 


for r=l, 2, ...,R. 

The equivalent charge centers for the case of R = 3 are computed by 
the following formulas as shown in figure 5: 


0,1 


bR 

V7 


0,2 


bR 


r - ^ 

a/5" 0,eq,2 


rQ 3 - bR 


(BlOb) 


0,eq,3 


bR 

/I 


In dealing with the beam loading problems, we require a complex ex- 
pression for the beam current, and furthermore, our main concern is the 
fundamental beam current. Thus, by letting n = 1 in the Fourier series 
expansion formula, we obtain 

ib(Zs’^s^ = Aq + Ai cos <(> + B^ sin (j> 



(Blla) 


with 


A|^ = arctan 
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and the fundamental beam current ibl(^S’^s) given by 


/ 2 2\^^^ 

- \''l * V '=OS(t - A^) 


= Re exp(j>) 


(Bllb) 


where 


Ibi = exp(-jAj^) (Bile) 

is the required complex beam current expression. 

Finally it is convenient to normalize the complex beam current to the 
dc beam current Ig. When this is done, we have the normalized beam cur- 
rent in complex form: 






-J^b 

e 


(B12a) 


where 




/2 


and 


(B12b) 


A 


b ~ 


arctan 



(B12c) 


With the beam current determined, the beam loading effect can now be 
determined by computing the beam admittance Yb as follows: 


Y 


b 






-j(e+Ajj) 

e 


(B13) 


where 
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complex beam current given by eg. (Bile) 

Vg complex gap voltage given by eq. (D6a) 

Gq dc beam conductance, Iq/Vq 

a gap voltage modulation index, |Vg|/Vo 

normalized Fourier coefficients, given by eq. (Bl2b) 

Ab phase factor of complex beam current, arctanCBi/Ai) , given by 

eq. (Bl2c) 

e phase factor of induced gap voltage, arctanC^l/U”!) - arctan(2Q-|-6j) , 

given by eq. (D6b) 

Finally the beam conductance G^, and beam susceptance can now be 
obtained by taking the real and imaginary part of Yb to yield 

^b = ^b = ^^1 ^1 ^b) 


and 


B 


b " 



(B15) 
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APPENDIX C 


DERIVATION OF EQUATIONS FOR CALCULATION OF INDUCED CURRENT 

The calculation of the induced current is based on the general formula- 
tion, equation (37), derived in a previous report (ref. 1). This is re- 
written in the following manner; 


where 

Pj(r,z,t) 

LTj(ro,zo;r,z,t) 



all j 


'cct, j 


TT 


ind 



(Cl) 


charge density of charge element 

vector velocity of charge element, which is 

function of entrance variables rQ and tQ, as 
well as its position (r,z) inside interaction gap 
at time t 


Ecct, i('"0’^0’'"’^’^) electric field induced in cavity by motion of charge 

element APj (=pj Atj) inside interaction region, 
which is defined (fig. 8) as linear (extended) gap 
distance from z$ - (1 + 2a) to Zg + (1 + 2a), 
where Zg is referred to midplane of cavity gap, 
s^h cavity 


'^ind 


induced voltage across cavity gap 


The induced voltage V-jp^ is related to the cavity induced electric field 
amplitude Eq by the following relation (see eq. (13)) of refs. 4 and 5): 


E 


H V 


ind,s 
(2 sinh Ha) 


(C2) 


where IVipd^gl is the induced gap voltage of the s'^h cavity, and H 
is the field-shape parameter. For a uniform field, H = 0. By the charge 
conservation principle, we can relate the charge element Aqj in any 



deformed disk charges at various positions inside the gap at time t ■ tj. 
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arbitrary interaction gap, the cavity gap, to the charge element 

AQj (=Pio Atj) at the entrance of the input cavity, where the beam is 
unmodulated and has a constant charge density Piq. This, in turn, can be 
expressed as follows: 


Q 0 
= - W = - W¥ = 

where Iq is the dc beam current, f is the frequency of the klystron, 
and N and R are the total number of disks and rings used in the computer 
model. With this transformation, equation (Cl) can be written as follows: 


NR 


’-TiW ^ Uj(r„,t„;r,z,t) 


. ‘-cct,j''0’''0 


(rn,tn;r,z,t) 


j=l 


where 


2 - {i + 2a) < z < z + (j, + 2a) 
s — — s 

This can be conveniently expressed by using Kronecker delta notation defined 

by 


6 




r 1 for z = 2 ^ 
\o for z ^ 2|^ 


where z^ is defined by 

- {i + 2a) 1 Z|^ < Zg + (*. + 2a) 


Thus 


^ ind 


(t) 




Uj(r,t;r,z,t) 


’cct, j 


(r,t,;r,z,t) 


(C3) 


all k j=l 

Expressed in terms of its Fourier components, equation (C3) can be expanded 
as 


'ind("k’^) = ^0 Z (^n °n 

n=l 
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and, in particular, the fundamental component of the induced current, 
■iind,l is given by 


'lnd,l - <='>*(♦ - ‘i) 


-JA, 


= Re 




= Re I . . , e 
ind,l 




(C5) 


where Iind,l is the complex expression for the induced current and is 
given by 


^ind 




(C6a) 


and 


-1 / 

^i = 


(C6b) 


in which Ci and Di are the Fourier coefficients given by 


■)- f 


2ir 


''ind,l^^k’'*^^ 


cos 


sin ♦ 


>0 


A(j)„ all k 
m 


NR 


Z Z Z 


j=l 




“ 


'cos <f 


m 


A4> 


sin 4> 


m 


m , 


36 



M=36 


NR r 


- w(r) X X '^-"k X 


m=l all z, 


j=l 


Uj(ro,t(,;r,Z|^t) ^ 


'cos <|| 


m 


(C6c) 


sin 4i 


m . 


In this equation, Ai|>ni has been replaced by (2 it/M) and M is taken to 
be 36. Usually the fundamental induced current is normalized in terms of 
the dc beam current Iq: then 


I, 


I, 


ind,l I 


i nd, 


0 


1 . 


H? ^ “^'^i 

i g I 


(C7a) 


where 


and 


. -1/ ^1 
A- = tan 
1 ' r 

^1 


(C7b) 


^i = t: 


(C7c) 


" T 


0 


(C7d) 


Finally the dot product inside the bracket can be expressed in terms of 
its two scalar components in the form 

E 

^ " V = T~ ^zY~ 

Using equation (C2), we obtain 


7 T . cct, j 

3 2 sinh (Hi) 
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By substituting this relation into equation (C6c), we get 


M=36 



2TrI^Ha 

inhflHa ^ ^ 


NR=(32)(3) 


D^/ MNR sinh 


(h=I) 


m=l all k 


j=l 




fcos 


(C8) 


where 


= 


r,3 

ao) 


and 


are the normalized velocity components in the radial and axial directions, 
respectively. They are functions of both the entrance position and phase as 
well as the position at the interaction gap and the time. 
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APPENDIX D 


DERIVATION OF EQUATIONS FOR DETERMINATION 
OF VOLTAGE MODULATION INDEX o 

With reference to figure 9 and reference 6, we see that the input 
admittance of an arbitrary cavity is given by 


" jS: " 


(Dl) 


where 


L cavity equivalent inductance 

C cavity equivalent capacitance 


G(. cavity conductance representing losses in a cold cavity, l/Rg 

G]_ equivalent load conductance as seen by cavity output window (perfect 
matched condition assumed), l/R^ 

Yt) beam admittance, Gb + 

Gb beam conductance. Re Yb 

Bb beam susceptance, Im Yb 


If we designate Gj = G^ + Gb + Gl as the total shunt conductance of the 
cavity, then, in the vicinity of resonance, when the cavity is tuned to the 
higher side of the resonant frequency fQ, by replacing lo = ioq + Aw, we 
can write equation (Dl) as follows: 

Y. = Gt + jB, + jw.C + j Aw C - j ^ (l + 

T D T D 0 ^ \ ^0 • 



+ j 



+ ^ Aw \ 

“o'- “0/ 


(D2a) 



Figure 9. - Generalized cavity equivalent circuit. 


39 



To obtain equation (D2a), we used the following relations: 


1 


+ -T — j- = 0 

0 Jo)qL 


WqLC = 1 


Furthermore, by denoting 

, Aio 


and 


R 





we can write equation (D2a) as 





1 + j2Qi 



+ 1 R R 

7 % U 


)] 


(D2b) 


where Qj is the total or overall cavity Q defined by the following 
relation : 


1 _ I r _ 1 + 

^ 



(D3) 


where 

Ql, Q of the cavity that would result if the cavity were cold (i.e., in 
the absence of the beam) and no external load were coupled to it, 
Rc/cogL, with or R^ representing losses in the cavity 

Qe Q that would result if the cavity were loss free and only loading by 
the external load were present, Rl/wqL 

Qb Q of the cavity due to beam loading alone, l/uoLGb = l/Gb(R/Q) 

Next let us introduce the new parameter 6y defined by 

^ ^ ^^drive “ ^res^ 
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where 6^J is the cavity frequency-detuning parameter in the absence of 
the beam, and (1/2)B|3(R/Q) is the cavity frequency-detuning parameter due 
to the contribution of the beam susceptance In terms of this new 

frequency-detuning parameter 6j, we write equation (D2b) as follows: 

= G^(l + j2Q^6^) 


= |/l + (2Q^6^)^ arctan ^(2Q^5-p) 


(D4b) 


With the input admittance of the cavity determined, the RF gap voltage 
can be found as follows: 


Vg(„) . 


^ i nd , 1 


(D5) 


where is the cavity induced current (fundamental) given by equa- 

tion (C6a),’and Y-j^(o)) is the input admittance of the cavity given by 
equation (D4b). When these relations are substituted into equation (05), we 
get the complex gap voltage 


r p 

V = Qt arctan e (D6a) 

g / 2* Q T 

yi + ( 2 q ^ 6 ^)^ 

where we have replaced Rj (= 1/Gj) by (R/Q)Qj, and the phase angle of 
the voltage is given by 


e 



- tan ^(2Qy6-|-) 


(D6b) 


Cj and are given by equation (C6c), and Qj and 67 are given by 

equation (03) and (04a), respectively. 

Finally the voltage modulation index a can be obtained by taking the 
modulus of Vg divided by the dc beam voltage Vg to yield 


a 



GqQt 


R 

Q 


+ 0^ 


where Gg = Ig/Vg is the dc beam conductance, and Ci and 0^ are 
the normalized Fourier coefficients given by equation (C7). 


(07) 
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